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ABSTRACT 


A circulation  control  uncambered  elliptic  airfoil 
section  with  a thickness-to-chord  ratio  of  0.20  was 
tested  subsonically  to  determine  its  aerodynamic 
characteristics.  Lift  coefficients  up  to  5 were  pro- 
duced at  momentum  coefficients  of  0.24.  The  initially 
high  unblown  drag  coefficients,  characteristic  of 
bluff  trailing  edge  airfoils,  were  greatly  reduced  at 
low  values  of  momentum  coefficient.  It  was  therefore 
possible  to  produce  equivalent  lift-to-drag  ratios 
in  excess  of  30  when  = 1.0.  The  ability  to  produce 
high  lift  coefficients  essentially  independent  of 
angle  of  attack  is  indicated  by  the  results  of  this 
investigation. 


INTRODUCTION 


Tangential  blowing  over  the  bluff  trailing  edge  of  a 20-percent 
uncambered  ellipse  was  investigated  experimentally.  This  configuration 
is  one  of  a continuing  series  of  two-dimensional  circulation  control 
airfoil  sections  ^ ^ which  employ  the  Coanda  effect  to  obtain  high  ii*lt 
augmentation  by  tangentially  ejecting  a sheet  of  air  near  the  trailing 
edge  on  the  upper  surface.  Because  of  the  Coanda  effect,  the  jet  sheet 
remains  attached  to  the  rounded  trailing  edge  and  provides  a mechanism 
for  boundary  la,er  control.  The  blowing  can  be  thought  of  as  a movement 
of  the  effective  aft  stagnation  point  thereby  producing  an  increase  in 
circulation.  For  further  discussion  of  this  circulation  control  concept 
and  its  potential  applications,  the  reader  is  referred  to  other  studies 

. ->  6"10 

in  the  series. 

The  results  of  this  investigation  are  being  used  to  ascertain  the 
effects  of  camber  on  the  aerodynamic  characteristics  of  a two-dimensional 
circulation  control  airfoil  section  by  comparing  the  data  in  this  report 
to  that  presented  earlier  by  Williams  and  Howe.-'-  Wind  tunnel  data  were 
generated  during  July  1974. 
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MODEL  AND  TEST  APPARATUS 


The  fiberglass  model  is  a modification  of  one  previously  tested  and 
is  based  on  an  analytical  ellipse  with  an  8-inch  major  axis  (chord)  and 
a 20-percent  thickness-to-chord  ratio.  The  rounding  of  the  trailing  edge 
to  produce  a 0.31- inch  radius  reduced  the  chord  to  7.908  inches.  An 
upper  surface  tangential  slot  consisting  of  a knife-edged  steel  blade  was 
placed  at  the  97-percent  chord  position.  The  slot  exit  was  the  throat 
of  a converging  nozzle  whose  height  was  adjustable  through  the  use  of 
pitch  screws.  An  undercut  was  made  on  the  blade  in  the  vicinity  of  the 
slot  to  ensure  that  flow  would  exit  tangentially  to  the  model  surface 
(see  Figure  1) . High  pressure  internal  plenum  chambers  and  flow  fences 
were  also  incorporated  into  the  model.  Tables  1 and  2 give  coordinates 
for  the  model. 

The  two-dimensional  tests  were  conducted  in  the  15  x 20  inch  subsonic 

tunnel  with  a vented  test  section  and  plexiglass  walls.  The  plexiglass 

permits  flow  visualization  be  means  of  tuft3  ard  oil  day-glo  pigment 

solution.  The  model  was  pressure  tapped  as  siiown  iu  Table  2.  Lift  and 

pitching  moment  coefficients  were  obtained  by  numerical  integration  of 

pressure  tap  readings  as  recorded  on  a multiple-por'.  scanivalve  readout 

system.  These  coefficients  were  corrected  by  the  addition  of  jet  reaction 
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components.  Standard  solid  blockage  corrections  were  applied  to  the 
measured  freestream  dynamic  pressure;  no  wake  blockage  factor  was  used 
because  of  the  uncertain  effects  of  the  jet. 

Drag  measurements  were  made  by  using  a drag  rake  placed  1.5  chord 

lengths  downstream  of  the  model  inclined  at  15  degrees  to  the  freestream. 

The  rake  employs  54  total  and  8 static  pressure  tubes,  with  the  heaviest 

concentration  of  tubes  near  the  center  height.  The  momentum  deficit 

12 

methods  of  Betz  and  Jones  were  then  used  to  determine  the  drag  coefficient. 
An  addition  of  mV^/q^S  was  then  made  to  the  drag  coefficient  to  account 
for  the  additional  jet  momentum. 

Past  practice  has  been  to  employ  separate  plenums  with  high  pressure 
air  sources  at  the  junction  of  the  model  and  tunnel  wall  to  prevent  wall 
boundary  layer  separation.  Although  the  present  model  was  equipped  with 
such  plenums,  it  was  decided  after  careful  analysis  of  prelir-iinary  data 
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that  they  would  not  be  used  because  of  damage  previously  sustained  by  the 
model.  A series  of  model  spanwise  taps  located  near  the  trailing  edge  was 
utilized  to  record  the  lateral  pressure  distribution  which  was  carefully 
monitored  during  the  test  run  as  an  indication  of  two-dimensionality. 

The  corrections  to  the  resulting  data  using  potential  flow  will  be  discussed 
later  in  the  report. 

Mass  flow  rate  (m)  was  measured  by  a calibrated  orifice  plate  inserted 
in  the  supply  line.  To  ensure  high  accuracy  for  low  mass  flow  data,  a 
second  orifice  plate  with  the  appropriate  calibration  range  was  used  for 
the  lower  values  of  C^.  The  jet  velocity  was  calculated  by  assuming 
isentropic  expansion  from  duct  stagnation  pressure  to  the  freestream 
static  pressure  as  follows: 


Vj  ajMj  = (YRTj 


)Js  M„  = 2 


2RT  (- 
t Y 
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The  momentum  coefficient  was  then  defined  as  = (mV^/qJS).  The  limiting 

value  of  C was  determined  by  consideration  of  the  maximum  safe  internal 
P 

plenum  pressure  for  this  particular  model  and/or  that  pressure  which 
would  result  in  jet  impingement  on  the  tunnel  floor.  Such  impingement 
could  cause  separation  of  the  floor  boundary  layer  as  well  as  inaccuracies 
in  the  drag  rake  data. 

A series  of  runs  was  made  at  freestream  dynamic  pressures  from  10  to 
30  psf  (Figure  2)  corresponding  to  a model  Reynolds  number  range  from 
2.4  x 105  to  4.2  x 105.  No  significant  effect  of  Re  on  the  data  was  noted 
and  = 20  psf  was  chosen  to  allow  for  a wider  range  of  C^. 

RESULTS  AND  DISCUSSION 


The  characteristics  of  the  20-percent  uncambered  ellipse  were  evaluated 
for  a single  slot  height  of  0.01  inch  (slot  height/chord  = 0.00126), 
momentum  coefficient  ranging  from  0 to  0.26,  and  angles  of  attack  a 
ranging  from  -20  to  8 degrees.  Figure  3 depicts  the  variation  of  momentum 
coefficient  with  duct  pressure  for  the  test  dynamic  pressure  of  20  psf. 
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The  expansion  of  the  slot  caused  by  the  pressurization  of  the  duct  is 
shown  in  Figure  4.  These  data  were  obtained  by  pressurizing  the  duct  and 
measuring  the  resulting  slot  height  with  a thickness  gage  under  quiescent 
tunnel  conditions. 


Figures  5 and  6 show  the  sectional  lift  coefficients  determined  by 
integration  of  the  midspan  taps.  Figure  5,  which  is  an  expanded  scale 
plot  including  closely  spaced  low  momentum  coefficient  data,  revealed 
the  possible  existence  of  a discontinuity  in  the  lift  slope.  The  density 
of  data  points  taken  in  previous  investigations  of  circulation  control 
two-dimensional  models  has  been  such  that  detection  of  a discontinuity 
in  the  lift  curve  slope  (dC./dC  ) would  be  unlikely.  However,  a change 
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in  the  slope  was  reported  by  Englar  . Oil  flow  studies  performed  during 
the  evaluation  of  this  model  indicated  that  the  flow  ahead  of  the  slot 
was  separated  under  conditions  of  no  blowing  (C^  =»  0) . This  discontinuity 
may  correspond  to  the  level  of  blowing  required  for  flow  attachment  to 
occur  at  the  slot  position.  Therefore,  the  existence  of  such  a discon- 
tinuity would  depend  on  the  chordwise  location  of  the  slot  in  relation 
to  the  flow  separation  point  with  no  blowing;  the  possible  dependency 
on  Reynolds  number  was  not  investigated  at  this  time. 

An  examination  of  Figure  6 shows  that  continually  increased  with 
increasing  C , indicating  that  "C^  stall"  had  not  yet  been  reached,  and 
it  appeared  that  the  lift  coefficient  could  have  been  increased  still 
further  with  increased  blowing. 

Another  characteristic  evident  from  Figure  6 was  a tendency  of  the 

lift  curve  to  "lie  over"  for  C >0.08,  resulting  in  a smaller  lift 

y 12 
coefficient  than  anticipated  in  view  of  results  reported  by  Englar  ’ . 

Additionally  a deviation  from  the  expected  lift  curve  was  experienced  at 

a = +8  and  +5  degrees.  This  phenomenon  was  particularly  pronounced  at 
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a = 8 degrees  and  has  been  noted  previously  ; however  in  this  case  a return 
to  predicted  values  was  never  obtained.  This  behavior  is  apparently  due 
to  the  aft  separation  experienced  on  this  model. 


Further  insight  may  be  gained  by  examining  figure  7 which  presents 

the  augmentation  ratio  as  a function  of  momentum  coefficient.  The  augu- 

mentacion  ratio  is  defined  as  AC^/C^,  where  is  the  increase  in  lift 

coefficient  above  the  unblown  value  for  a given  value  of  C and  incidence. 

y 

A significant  loss  of  augmentation  is  apparent  at  a ■ +8  degrees  and  to 
a lesser  extent  at  a = +5  degrees;  this  affirms  the  deviation  from  the 
expected  lift  curve  slope  noted  previously  in  Figure  6.  In  general  the 
augmentation  at  higher  values  of  momentum  coefficient  fall  off  faster 
than  with  previous  models,  it  should  be  noted  at  this  point  that  data 
for  a = 20  degrees  were  deliberately  omitted;  the  recording  of  accurate 
unblown  data  was  not  possible  and  therefore  a computation  of  augmentation 
ratio  could  not  be  made. 

With  few  minor  exceptions,  an  examination  of  plots  revealed  no 
evidence  of  leading  edge  separation  bubbles  although  minor,  local  upper 
and  lower  surface  trailing  edge  separations  were  common.  At  large  negative 
angles  of  attack  (a  * -20  degrees),  the  lower  surface  was  initially 
completely  separated,  but  reattachment  occurred  with  increased  blowing 
due  to  the  movement  of  the  leading  edge  stagnation  point  to  the  lower 
surface.  As  the  jet  velocity  was  increased  still  further,  the  lower 
surface  again  showed  areas  of  separation  in  the  vicinity  of  the  trailing 
edge. 

The  variation  of  lift  coefficient  with  geometric  angle  of  attack  is 
presented  in  Figure  8.  The  slopes  of  the  curves  were  similar  for  unstalled 
flow  conditions.  Agreement  between  the  theoretical  lift  curve  and 
experimental  results  was  good  for  conditions  of  no  blowing. 

To  complete  the  discussion  of  the  lift  characteristics  of  the  20-percent 
ellipse,  consideration  of  effects  induced  by  spanwise  nonuniformity  is 
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required.  The  suggested  use  of  flow  fences  and  tip  jets  (see  Englar  and 

8 

Englar  and  Williams  ) was  not  possible  in  this  case  for  reasons  previously 
discussed,  and  therefore  a determination  of  the  effective  angle  of  incidence 
was  made.  For  the  experimental  cases  selected,  potential  flow  pressure 
distributions  for  several  incidences  and  an  adjusted  were  produced. 

The  adjustment  to  the  lift  coefficient  required  that  the  increment  of 
lift  due  to  the  jet  suction  peak  be  determined  and  subtracted  from  the 
experimental  lift  coefficient.  Since  this  increment  would  not  be 
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theoretically  predicted,  these  distributions  were  then  compared  to  the 
experimental  pressure  distribution  until  the  leading  edge  stagnation 
point  coincided.  A resulting  experimental-potential  flow  match  is  shown 
in  Figure  9.  The  agreement  was  good  except  for  the  additional  jet  suction 
peak  and  minor  lower  surface  separation  (both  not  predictable  with  the 
use  of  potential  flow).  The  effective  angle  of  incidence  for  the  experi- 
mental data  is  depicted  in  Figure  10. 

The  value  oi  the  minimum  pressure  on  the  airfoil  as  a function  of 
lift  coefficient  is  presented  in  Figures  11  and  12.  These  data  are  mainly 
intended  for  marine  applications  of  circulation  control  systems,  where 
cavitation  is  a fun.';ion  of  minimum  pressure.  Figure  12  is  the  result 
of  the  low-blowing,  small  alpha  test  series  which  would  be  most 
appropriate  for  the  above-mentioned  application. 


Figures  13  and  14  indicate  the  variation  of  a modified  drag  coefficient 

with  momentum  coefficient.  These  data  result  from  an  integration  of  the 
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wake  deficit  using  the  methods  of  Jones  and  Betz  which  were  then  modified 
to  account  for  the  additional  momentum  of  the  jet,  thereby  becoming 
Cd  = cdrake  “ (oV^/qS).  Both  figures  indicate  a high  initial  drag  at  all 
alphas  primarily  because  of  the  nature  of  bluff  trailing  edge  airfoils. 

The  highest  initial  drag  was  at  a = -20  degrees;  visually  monitoring  by 
manometer  tubes  indicated  that  the  wake  was  extending  over  the  entire 
rake  before  a no-blowing  condition  could  be  reached,  and  this  was  attributed 
to  extensive  separation  on  the  lower  surface.  The  secondary  rise  in  the 
drag  curve  at  both  a = -15  and  -20  degrees  at  < 0.2  was  again  the 
result  of  lower  aft  separation  after  reattachment  had  occurred. 

The  initial  variation  of  drag  with  trailing  edge  blowing  can  be  more 
easily  observed  in  Figure  13.  An  initial  decrease  in  drag  was  immediately 
noted  whenever  a £ 0 an  initially  flat  drag  distribution  was  evident  for 
a = +8  or  +5  degrees.  This  was  also  evident  to  a lesser  extent  at  a = +3 
degrees  where  the  initial  decrease  in  drag  was  less  than  at  any  nonpositive 
alpha.  For  a = 8 and  5 degrees,  this  region  was  found  to  correspond  to 
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the  deficit  in  the  lift  curve  slope  previously  discussed.  Significant 
thrust  recovery  was  noted  in  Figure  14  for  most  alphas,  as  evidenced  by 
the  decreasing  drag  coefficients. 


I 1 


PITCHING  MOMENT 

Figure  15  depicts  pitching  moment  about  the  midchord  (Cm^) • As  was 
the  case  with  previous  circulation  control  sections,  a negative  pitching 
moment  was  produced  by  the  relatively  high  trailing  edge  suction  peaks. 

The  only  notable  exceptions  were  at  a = -15  and  -20  degrees  for  _>  0.134. 
In  both  those  cases,  there  were  large  regions  of  separation  on  the  lower 
surface  of  the  trailing  edge.  The  resulting  region  of  low  pressure 
reduced  the  negative  pitching  moment  somewhat. 


EQUIVALENT  LIFT-DRAG  RATIO 

Comparison  of  data  taken  with  and  without  blowing  can  best  be  made 
when  account  is  taken  of  the  energy  expended  to  produce  blowing.  Figure  16 
presents  the  lift-to-equivalent-drag  ratio  as  a function  of  lift  coefficient; 
the  equivalent  drag  is  defined  here  as: 

mV.2 

d = d + r-1-  + mV 
e 2V  00 


The  first  term  d is  the  momentum  deficit  as  measured  by  the  drag  rake 
(corrected  for  jet  efflux),  the  second  term  is  the  kinetic  energy  flux, 
and  the  third  is  an  intake  momentum  flux.  Thus  in  coefficient  form: 
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Cd  + Cp  2V 


+ C tt- 
P Vj 


A more  comprehensive  derivation  of  the  above  may  be  found  in  Englar  . 

The  highest  value  for  &/de  occurred  at  very  low  blowing  and  the  maximum 
value  was  approximately  40,  despite  the  relatively  high  lift  coefficients. 
Since  the  measured  value  of  drag  was  low  throughout  a great  portion  of  the 
experimental  test  range,  the  results  show  the  importance  of  "blowing  power" 


(C  V./2V  ) in  the  expression  for  equivalent  drag. 
V j 00 


Thus  to  obtain  high 


values  for  l/d  it  becomes  necessary  to  produce  high  lift  coefficients 
e 

at  relatively  low  values  of  the  blowing  coefficients. 


CONCLUSIONS 


The  20-percent  thick  circulation  control  ellipse  was  evaluated  in  a 
subsonic  wind  tunnel  to  provide  continuing  information  on  a new  family 
of  airfoils.  The  following  conclusions  were  drawn  from  an  examination 
of  the  subsonic  data: 

1.  Despite  the  fact  that  augmentation  ratios  were  generally  lower 
than  obtained  in  previous  test  of  circulation  control  airfoils,  section 
lift  coefficients  cf  5.0  were  generated  for  C^  = 0.24  at  an  angle  of 
attack  of  3 degrees.  A lift  coefficient  of  2 was  generated  by  the  section 
at  a = -20  degrees. 

2.  Augmentation  ratios  in  excess  of  75  (see  Figure  7)  were  produced 
at  low  blowing  rates  for  negative  angles  of  incidence. 

3.  At  C = 1.0,  lift-to-equivalent-drag  ratios  of  over  30  were 
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generated. 

4.  In  most  cases  initiation  of  blowing  reduced  the  high  unblown  drag 

coefficient.  In  those  cases  where  the  value  of  the  drag  coefficient 

remained  unchanged  until  a substantial  level  of  blowing  was  reached,  a 

low  initial  lift  slope  (dC./dC  ) was  also  indicated. 

jc  y 

5.  The  presence  of  a possible  discontinuity  in  the  ?.ift  curve  slope 

(dC  /dC  ) may  correspond  to  the  movement  of  the  separation  point  to  a 

36  y 

position  aft  of  the  slot  location. 
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TABLE  1 

Two-Dimensional  Model  Coordinates  for  Upper  and  Lower  Surfaces 


Table  1A 

- Upper  Surface 

Table  IB 

- Lower  Surface 

X 

Y 

X 

Y 

0.0000 

0.0003  LEADING  EDGE 

7.8980 

-0.0364  TRAILING  EDGE 

0.0200 

0.0682 

7.8780 

-0.0742 

0.0400 

0.1095 

7.8580 

-0.1010 

: „ 

0.0600 

0.1402 

7.8280 

-0.1340 

0.0799 

0.1648 

7.7980 

-0.1628 

0.1000 

0.1859 

7.7680 

-0.1861 

0.1300 

0.2103 

7.7380 

-0.2081 

0.1600 

0.2288 

7.6980 

-0.2339 

0.1800 

0.2429 

7.6580 

-0.2569 

0.2200 

0.2527 

7.6180 

-0.2774 

j 

0.2500 

o.2a*8 

7.5781 

-0.2972 

| 

0.2950 

0.2993 

7.5380 

-0.3152 

0.3400 

0.3140 

7.4930 

-0.3340 

0.3850 

0.3288 

7.4230 

-0.3611 

0.4300 

0.3429 

7.3630 

-0.3819 

! k 

0.4750 

0.3569 

7.3030 

-0.4017 

0.5450 

0.3784 

7.2530 

-0.4166 

0.7050 

0.4245 

7.0930 

-0.4595 

i 

0.8650 

0.4651 

6.9330 

-0.4963 

1.0250 

0.5010 

6.7730 

-0.5291 

1.2250 

0.5401 

6.6130 

-0.5584 

1.3850 

0.5672 

6.4530 

-0.5829 

l 

1.6050 

0.6008 

6.2930 

-0.6055 

2.0050 

0.6509 

5.8779 

-0.6551 

2.4050 

0.6893 

5.3780 

-0.7010 

2.8050 

0.7181 

4.8780 

-0.7328 

3.2050 

0.7367 

4.3780 

-0.7499 

3.6300 

0.7485 

3.8780 

-0.7545 

4.1300 

0.7522 

3.4030 

-0.7515 

4.5050 

0.7464 

3.0030 

-0.7403 

5.0051 

0.7257 

2.6030 

-0.7169 

5.5050 

0.6930 

2.2030 

-0.6828 

6.0050 

0.6428 

1.8030 

-0.6363 

6.3350 

0.6000 

1.4630 

-0.5851 

6.4950 

0.5771 

1.3030 

-0.5573 

* 

6.6550 

0.5520 

1.1430 

-0.5279 

* ' 

6.8150 

0.5238 

0.9830 

-0.4947 

6.9750 

0.4922 

0.8230 

-0.4569 

7.1350 

0.4557 

0.6630 

-0.4130 

7.2700 

0.4191 

0.5030 

-0.3622 

- 

7.3300 

0.4011 

0.4430 

-0.3419 

7.3900 

0.3814 

0.3830 

-0.3208 

5 " 

7.4500 

0.3595 

0.3229 

-0.2993 

| 

7.5100 

0.3357 

0.2630 

-0.2761 

* 1 

7.5500 

0.3180 

0.2180 

-0.2565 

7.5900 

0.2988 

0.1680 

-0.2311 

\ 

7.6300 

0.2779 

0.1380 

-0.2121 

\ 

7.6700 

0.2546 

0.1080 

-0.1898 

7.7100 

0.9292 

0.0780 

-0.1617 

7.7500 

0.2015 

0.0480 

-0.1236 

* * 

7.7900 

0.1696 

0.0280 

-0.0883 

1 

7.8201 

0.1418 

0.0080 

-0.0346  LEADING  EDGE 

* < 

7.8500 

0.1099 

7.8800 

0.0696 

3 

7.9000 

0.0281 

\ 

7.9080 

0.0000  TRAILING  EDGE 

l 
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TABLE  2 

Two-Dimensional  Model  Pressure  Tap  Coordinate  Listing 


' Surface 


Lower  Surface 


p No. 

X Position 

Tap  No. 

X Position 

1 

0 

55 

0.125 

2 

0.084 

54 

0.205 

3 

0.177 

53 

0.301 

4 

0.285 

52 

0.396 

5 

0.410 

51 

0.606 

6 

0.590 

50 

0.801 

7 

0.782 

49 

1.188 

8 

1.186 

48 

1.601 

9 

1.591 

47 

2.011 

10 

2.001 

46 

2.421 

11 

2.403 

45 

2.801 

12 

2.795 

44 

3.210 

13 

3.184 

43 

3.624 

14 

3.618 

42 

4.293 

15 

4.295 

41 

4.743 

16 

4.747 

40 

5.194 

17 

5.201 

39 

5.601 

18 

5.603 

38 

5.991 

19 

5.996 

37 

6.400 

20 

6.394 

36 

6.806 

21 

6.799 

35 

7.191 

22 

7.201 

34 

7.437 

23 

7.666 

33 

7.649 

24 

7.710 

32 

7.737 

25 

7.723 

31 

7.848 

26 

7.801 

30 

7.888 

27 

7.843 

29 

7.906 

28 

7.885 

29 

7.906 

M 
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